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I Abstract. In this paper we study a one-dimensional quantum random walk with the Hadamard 
transformation which is often called the Hadamard walk. We construct the Hadamard walk using 
a transition matrix on probability amplitude and give some results on symmetry of probability 
', distributions for the Hadamard walk. 

1 Introduction 

X. 

^ , Recently, quantum versions of classical random walks have been investigated by a number of 
groups, for examples, see refs. [1-10], [13-22], [24], [26-29]. A more general setting including 
quantum cellular automata can be found in ref . [23] . It is well known that classical random 
walks have found practical applications in many fields including physics, chemistry, biology, 
astronomy, and computer science. The hope is that quantum random walks might lead to 
applications unavailable classically. In fact, quantum random walks have many suggested 
applications in physics and quantum information. For examples, quantum walks have the 
potential to give various new tools for quantum algorithm design [3]. It is also believed 
that quantum walks can provide a benchmarking protocol for ion trap quantum computers 
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[28]. Childs et al. [9] shows query complexity of some problems is speed up exponentially by 
quantum walk. Kempe [13] and Yamasaki et al. [29] are the first which showed exponential 
speed up of hitting time in discrete time walk, while Childs et al. [9] treats continuous time 
model. (Yamasaki et al. [29] conjectured the speed up by numerical computation, and the 
conjecture was proven by Kempe [13].) So to clarify some properties of the distribution for 
the quantum random walk can be considered as important work. 

In the present paper, we focus on the one-dimensional Hadamard walk whose time evo- 
lution is given by the Hadamard transformation (see p. 19 in ref. [25]): 



V2 



1 1 
1 -1 



The Hadamard matrix H is unitary. The Hadamard walk considered here is a quantum 
generalization of a classical symmetric random walk in one dimension with an additional 
degree of freedom called the chirality. The chirality takes values left and right, and means 
the direction of the motion of the particle. The evolution of the Hadamard walk is given by 
the following way. At each time step, if the particle has the left chirality, it moves one step 
to the left, and if it has the right chirality, it moves one step to the right. 

More precisely, the Hadamard matrix H acts on two chirahty states \L) and \R): 



\R) 



^{\L) + \R)) 



V2 



i\L) - \R)) 



where L and R refer to the right and left chirality state respectively. In fact, define 



\L) 



1 






1 



so we have 



1 
1 



H\L) 

H\R) = ^(|L) - (it:)) 



Here we introduce P and Q matrices as follows: 



V2 



1 1 






1 -1 



with 



H = P + Q 



Here P (resp. Q) represents that the particle moves to the left (rcsp. right) with equal 
probability. We should remark that P and Q are useful tools in the study of the iterates of 
H. However, they cannot be interpreted as dynamical evolution operators since they are not 
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unitary. By using P and Q, wc define the dynamics of tlie Hadamard walk in one dimension. 
To do so, we introduce tlie next (2A^ + 1) x (27V + 1) matrix : {C'^f^+^ (c2)2iv+i . 

P Q 

Q P 

Q P ... 



H 



N 



... g p 

Q P 

P Q 



wliere C is tlie set of complex numbers and 

= 

Let 








(n) 



be the two component vector of amplitudes of the particle being at site k and at time n with 
the chirality being left (upper component) and right (lower component), and 

*W(^) = • • • , ^Sv"^^^)] e (c^)^^+^ 

be the qubit states at time n where t means transposed operator. Here the initial qubit state 
is given by 



N 



N 



*(°)(^) = *[C^, v:',C^ e (C^)^^^^ 



where 







" " 




a 










with lap + = 1 and is the absolute value of z e C. The following equation defines 
the time evolution of the Hadamard walk: 

where (*(")((/?))fe = ^k'\v)- Note that P and Q satisfy 

PP* + QQ* = P*P + Q*Q = /, PQ* = QP* = Q*P = P*Q - 








where * means the adjoint operator. The above relations imply that becomes also 
unitary matrix. 



3 



N N 

For initial state Tp — *[oy~?"^, Lp,'0^^^^~^, we have 

N-l N-l 

HnTp = *[or^, P^, 0, g<^, 



Ar-2 



Ar-2 



Hlif = *[0r^, PV, 0, {PQ + QP)<f, 0, QV, 



Ar-3 



Ar-3 



i/^r^ = *[0, . . . , 0, PV, 0, (P'Q + PQP + QP')^, 0, (Q'P + QPQ + PQ')ip, 0, QV, 0, . . . , 0] 

This shows that expansion of if" = (P + Q)" for the Hadamard walk corresponds to that of 
2" = (1 + 1)" for symmetric classical random walk, that is, Pascal's triangle. 
Define the set of initial qubit states as follows: 





a 


1- 





e C 



2- 'aP 



Now we define probability distribution of the Hadamard walk X!^ starting from initial qubit 
state </? e $ by 

p«=fc) = i*i"V)r 

The purpose of this paper is to study some properties of distribution of X!^, for example, 
the symmetry of distribution. To do so, we introduce the following classes of initial qubit 
states: 



$0 





a 


1- 





\a\ 



I, al3 + a/3 = 



{(pe^: P{X^ = A;) = P{X^ = -k) for any n e Z+ and A; e Z} 
e $ : E{X:^) = for any n e Z+} 



and Z (resp. Z+) is the set of (resp. non-negative) integers. We should remark that if 
q;/3 7^ 0, then a/? + a/3 = imphes that there is a A; e Z such that 9a — Op + (7r/2) + kn, 
where z is conjugate of z G C and 6z is the argument of z G C, that is, a and /3 are 
orthogonal. For ip G the distribution of is symmetric for any n G Z_|_. 
In this setting, we have the next result (Theorem 2 in Section 2): 



This result can be generalized to a general quantum random walk in one dimension (see 
Konno [17, 18], for example). However here we give other two proofs which appear in 
Sections 2 and 3 respectively. 

The rest of the paper is organized as follows. In Section 2, we study the necessary and 
sufficient conditions of the symmetry of distribution of X!^ for any time step n. Section 
3 treats another approach to the study of the symmetry by using a quantum version of 
Pascal's triangle (see Theorem 3). Section 4 is devoted to a conjecture on the expectation of 
the Hadamard walk. Appendix A gives a proof of Theorem 3. In Appendix B, we compute 
the mth moment of a limit distribution for the rescaled Haramard walk in the symmetric 
case. 
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2 Symmetry of Distribution 

First we present the following useful lemma to prove Theorem 2. Let 



J 



-1 

1 



Then we have 

LEMMA 1 (i) We suppose that initial qubit state is 



V2 



where 6 e [0, 27r) . Then we have 

(2.1) - i-mj'^'^-li^p) 

for any k & Z and n e Z_|_ . 

(a) We suppose that initial qubit state is 



1 

—i 



where e [0, 27r) . Then we have 



(2.2) 

for any k & Z and n e 



in), 



Proof. We show parts (i) and (ii) by induction on time step n. The proof of part (ii) is 
essentially the same as that of part (i), so we omit it. For simplicity, we write \E'^"''* = '^^^\lp). 

First we consider n = case. For any k ^ 0, the initial condition gives = 0, so we 
can check Eq. (2.1). As for /c = 0, we have 



(0) 



" -1 " 




' 1 ' 


1 


72 


i 



(0) 



Therefore Eq. (2.1) is correct for n — 0. 

Next we assume that Eq. (2.1) holds for time n — m. We should remark that the 
definition of the Hadamard walk gives 



(2.3) 

for any k. Moreover note that 



(2.4) 



QJ = - JP, PJ = -JQ 
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Then we have 



(m+l) 



= g(-i)"^ij^ 



(m.) 
(fc-1) 



(m) 

-(fc+1) 



(m) 



(m+l) 



The first and fifth equalities are given by Eq. (2.3). The second equality comes from the 
induction hypothesis. The fourth equality is obtained by Eq. (2.4). So it is shown that Eq. 
(2.1) is also correct for time n — m + l. The proof of Lemma 1 is complete. 

The following three classes were given in Introduction. 





a 


1- 





e $ : b| 



a/3 + a/3 = 



= {(^ e $ : |*i"H(^)| = |*L"^((^)| for any n e Z+ and A; e z} 
$0 = {cpe^: E{X'^) = for any n e Z+} 

As we mentioned in Introduction, it is noted that a/? + a/3 = with a(5 ^ implies that a 
and /3 are orthogonal. For (/? e the probability distribution of is symmetric for any 
n e Z+. By using Lemma 1, we obtain the following result: 

THEOREM 2 For the Hadamard walk, we have 



Proof. First we see that the definitions of $s and $o give 
(2.5) C $0 

Next we should remark that 



$1 





" 1 ' 




1 


172 


i 


' 72 


—i 



:^i,^2e[0,27r)L 



By the above remark and Lemma 1, for any e we have 

\^t\^)? = (-i)^"Kr*L1(^)v*M1(^) = 

Therefore if (/? e then 

for any n e Z and /c e Z+, so we have 

(2.6) $x C 
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A direct computation gives 



E{X^) ^ E{X2) = -{aP + aP) 

The above equations imply that if </? e $o then (/? e So we have 
(2.7) $0 C $± 

Combining Eqs. (2. 5)- (2. 7) gives 

$± C $s C $0 C $± 
so the proof of Theorem 2 is complete. 



3 Quantum Pascal's Triangle 

In this section, we take another approach to the study on symmetry of distribution of X!^, 
more precisely, distribution of P{X!^ — k) ior n + k — even. To do so, we need to consider 
the following expression. For fixed I and m with I -\-m — n and m — I — k, 

E{1, m) = Yl P^^Q'^'P^^Q"^^ ■ ■ ■ P^-g*"" 

Ij ,mj >0:mi H \-mn=Tn,h-\ \-ln=l 

For example, in the case of P{X4 = —2), we need to know the expression: 

S(3, 1) = P^Q + P'^QP + PQP"" + QP^ 
We remark the next useful relation: 

By using these, S(3, 1) becomes 

s(3, 1) = ^pg + -^PQP + ^PQP + \qp 



Here we introduce other useful matrices: 



R^V2PQ 



V2 



1 -1 




, S = V2QP 



V2 





1 1 



In general, we obtain the next table of computations with P, Q, R and S: 



7 





P 


Q 


R 


s 


p 


P 


R 


R 


p 


Q 


s 


-Q 


Q 


-s 


R 


P 


-R 


R 


-p 


S 


s 


Q 


Q 


s 



where we omit l/\/2 factor, for example, PQ = R/\/2. Since P, Q, R and S are a basis for 
the set of 2 X 2 matrices with complex valued components, we have 

E;(/, m) = pn{l, m)P + m)Q + r„(Z, m)R + m)S 

Next problem is to obtain explicit forms of Pn{l,rri),qn{l,rn),rn{l,'m) and Sn{l,m). In the 
above example oin — l + m — A case, we have 

S(4,0)= ('-^'j P, E{3,1)^(-^) {2P + R + S\ S(2,2)= f-^j (-P + g), 



S(l,3)=(^i=j {-2Q + R + S), S(2,2) = -(^i=j Q 



For a general S(Z, m), the following result is obtained. 



THEOREM 3 We consider the Hadamard walk. Suppose that l,m > with I + m — n, 

then we have 

(i) for I A m{— min{Z, m}) > 1, 
(a) for l{— n) > l,m = 0, 

H(/.0)=(i^)"p 

(Hi) for I = 0, m(= n) > 1, 

2(0,m)= (^-^) {-ir^'Q 

Theorem 3 can be obtained as a special case of Lemma 1 in Konno [18]. However for the 
convenience of readers, we give the proof in Appendix A. 

Prom now on we will give an application of this theorem as follows. By Theorem 3, we 
can get C For Z A m > 1, 

*S(Z,m)S(Z,m) 



7 



P- 



m 



7 



7 



7 



Q + R + S 
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D""'i:gM)-i;:l)(::l)(::.')(::.^: 



7=1 5=1 

'l-S 



P-Hl^Q + R + S 
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7 



n lAm lAm 

EE 



.l)7+'5 // _ 1\ // _ 1\ /„j _ A 1^ 



X 



" ;2 + (^_27)(m-25) /(/- 27) - m(m - 25) 
- 2(5) - m(m - 27) + - 27)(/ - 2(5) 

To consider the symmetry, for Z A m > 1, we have to check 

m)E{l, m) - ''E{m, l)E{m, I) 

2) hh~^\i-i)\5-i)U-i)\5-i^ 



7=1 5=1 

x2(Z -m) 

Similarly, for / > 1 and m = 0, we have 



7 + (5 n-2(7 + 5) 
n-2(7 + (5) -(7 + (^) 



n— 1 





n-1 


" 


1 " 


{I 


) 


1 






*5(/,o)s(;,o)-*s(o,os(o,/) = 

In any case, for any Z, m > except I — m — there exist a, 6 e C such that 

*5(i, m)H(Z, m) - *5(m, Z)S(m, Z) 

Therefore we obtain 



a b 
b —a 



in), 



a b 




a 


b —a 







a{\a\'^ - \P\'^) + b{aP + aP) 



where k = m — I. Prom the above result, we see that if (/? = /3] G then (/? e $s 
immediately. 



4 Discussions 

In the case of the Hadamard walk, a direct computation implies that for n G {1, 2, . . . , 10}, 

= -a„(|a|2 - |/?|2) - bn{a^ + aP) 

where 

1 9 5 27 17 293 157 

oi — 02 — 0, 03 — -, 04 — 1, 05 — g' "6 — ^, 07 ~ Yg' ~ "g"' ^9 ~ Y28' ~ "eT 

3 , ^ , 17 , 9 , 43 , 25 , 421 
61 = 62 = 03 = 1, 04 = 2' ^ ^ ' ^ ^ ^ ^ 4' ^ ^16' "s"' " 128 



So we conjecture that bn+i = a„ + 1 for any n > 1. This property would be useful for 
studying symmetry of distribution. 

Recently Konno [17, 18] gave a new type of limit theorems for quantum random walks 
in one dimension. In the symmetric case, his result implies that X!^/n converges weakly to 
Z"^ whose mth moment (= c^) is given by 

E{{Z^f-) = 1 - ^ E ^ (7 j > E{{Z^r-^) = 

for any n > 1 (a derivation of the mth moment can be seen in Appendix B). By using finite 
moments (m = 0, 1, . . .), we obtain the associated orthogonal polynomials {Pm}. So the 
interacting Fock space corresponding to the limit distribution for our symmetric Hadamard 
walk can be derived from the usual way (for example, see pp. 73-76 in Hashimoto, Hora and 
Obata [12]). To clarify the relation between the quantum random walk and the quantum 
decomposition is one of the future interesting problems. 

Appendix A: Proof of Theorem 3 
Here we give a proof of Theorem 3. 

(a) pn{l,m) case : First we assume I > 2 and m > 1. In order to compute Pn{l,m), it is 
enough to consider only the following case: 

iri W2 wy, 'W2-i W27 + I 

C{p, {I, m) = PP^QQ^~QPP^ ■ ■ ■QQ^~QPP^^ 

where w — (wi, W2, ■ ■ ■ , W2-y+i) G Z'^~^^ with wi, W2, • • • , W2.y+i > 1 and 7 > 1. For example, 
PQP case is Wi — W2 — W3 — 1 and 7 = 1. We should remark that I is the number of P's 
and m is the number of Q's, so we have 

I ^ Wi + Ws -\ h 

m = W2 + W4-\ h W27 

Moreover 27 + 1 means the number of clusters of P's and Q's. Next we consider the range 
of 7. The minimum is 7 = 1, that is, 3 clusters. This case is 

P---PQ---QP---P 

The maximum is 7 = (/ — 1) A m. This case is 

PQPQPQ ■ ■ ■ PQPQP ■P (Z - 1 > m), PQPQPQ ■ ■ ■ PQPQ ■■■QP {l-l<m) 

for examples. Here we introduce a set of sequences with 27 + 1 components: for fixed 

7 e [1, (/ - 1) Am] 

W{p,2j+1) = {w ^ {wi,W2,- ■ ■ ,W2^+l) ^ ■■ Wi+W3-\ hW27+l = ^, 

ty2 + W4 H h W2-y = m, wi, W2,..., W2j, ^27+1 > 1} 
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From the next relations: 



-P^ — — '7=Pi — — '7=Qi PQ — — 7=R-> — — 7=Rj RP — — i=P 

^2 ^2 V2 V2 V2 



we have 

where ty e 27 + 1). For Z > 2, and m > 1, that is, 7 > 1, we obtain 

Note that the right hand side of the above equation does not depend on w G W{p, 27 + 1). 
Finally we compute the number of w = {wi, W2, ■ ■ ■ , W2-y, W2-y+i) satisfying w £ W{p, 2j + 1) 
by a standard combinatorial argument as follows: 

|W^(P,27+1)|= ^^"^^ ^"""^^ 



7 / Vt- 1. 



From the above observation, we obtain 



(/-l)Am 

Pn{l,m)P = E E Cip,w)^:'y^'\l,m) 

7=1 «;6W^(p,27+l) 
(Z-l)Am 

= E \W{p,2j + l)\C{p,wY^^^'\l,m) 



7=1 \ 

Therefore the desired result is obtained: 



1 \ , . „, ntA . — l\/m — 1 
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When I > 1 and m — 0, this case is 

\p - pi - ( 



SO we have 

P»((.0)=(^)" 

Furthermore, when / = 1, m > 1 and / = 0, m > 0, we see 

Pn{l,m) = 

(b) qn{l,m) case : As in the case of part (a), we assume I > 1 and m > 2. In order to 
compute qn{l,m), it is sufficient to consider only the following case: 

Ull W2 UI3 W2~f+1 

C{q, wY^^+'^ {l,m)^ QQ^~Q'pf^QQ^~Q ■ ■ ■ Pp'^QQ^~Q 

where w — (wi, W2, ■ ■ ■ , ^27+1) £ Z'^'^^ with wi,W2,..., ^27+1 > 1 and 7 > 1. For example, 
QPQ case is = ^2 = ^3 = 1 and 7 = 1. Then we have 

I ^ W2 + W4 -\ h W27 

m^Wi + W3-\ h W27+1 

As in the previous part, the range of 7 is between 1 and / A (m — 1). Here we introduce a 
set of sequences with 27 + 1 components: for fixed 76 [1,1 A {m — 1)] 

W{q,2'y + 1) = {w = {wi,W2,- ■ ■ ,W2^+i) e Z^^^^ : W2 + W4^-\ \-W2^ = l, 

Wi+W3-\ h W27+I = m, Wi, ^2, . . . , -^27, -^27+1 > 1} 

Similarly we have 

for Z > 1, m > 2 and w e 27 + 1). Furthermore 
From these, we obtain 



\W{q,2j+l) 



{l-l)Am 

qn{l,m)Q = E E C{q,wY^^^'\l,m) 

7=1 «,eW(q,27+l) 

;A(m-l) 

= E |W(g,27+l)|C(g,^/;)(f^+^)(/,m) 

= (S)'-.-.,--f".-...(;:;)(-;> 
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Therefore the desired result is obtained: 

...,..,=-(^)'-V,r-|",-.,.(;;;)(",-) 

When / = and m > 1, this case is 

qn{0,m) = 

Furthermore, when Z > 1, m = 1 and Z > 0, m = 0, we see 

qn{l,m) = 

(c) rn{l,m) case : First we assume I > 1 and m > 1. In order to compute rn{l,m), we 
consider only the following case: 

Wl W2 W3 W2-/ 

C{r, «;) (f^) {I, m) = 'PP^~PQQ^TP^~P ■ ■ ■ QQ^~Q 

where w = {wi, W2, ■ ■ ■ , W27) G with wi, W2, ■ ■ ■ , W27 > 1 and 7 > 1. For example, PQPQ 
case is Wl — W2 — W3 — W4 — 1 and 7 = 2. Note that 

/ = Wl + W3 + h 

m = W2 + W4 H h 

Moreover 27 means the number of clusters of P's and Q's. Next we consider the range of 7. 
The minimum is 7 = 1, that is, 2 clusters. This case is 

P--PQ--Q 

The maximum is 7 = Z A m. This case is 

PQPQPQ---PQP---PQ (l>m), PQPQPQ ■ ■ ■ PQPQ ■ ■ -Q (m>l) 
for examples. Here we introduce a set of sequences with 27 components: for fixed 7 e [1, 1 Am] 
l^(r, 27) = {w ^ (wi,W2,- ■ ■ ,^2^) e Z'^'^ : wi + W3-{ \-W2-y-i^l, 

W2 + W4-\ h W27 = m, Wl, W2,..., W2J-1, W27 > 1} 

As in the case of part (a), Eq. (3.9) becomes 

for Z > 1, m > 1, that is, 7 > 1, and w e W{r, 27). Furthermore 



\W{r,2^)\^ 





From the above results, we obtain 



l/\m 

rn{l,m)R = E E C{r,wY^^\l,m) 

7=1 weW{r,2'y) 
I Am 

- J2\W{r,2j)\C{r,wY^''\l,m) 

7=1 

n-1 



Therefore 



.(.™)=a)""(-r|(-r(;:;)(::0 

When lAm — 0, this case is 

r„(/,m) = 

(d) s„(/,m) case : First we assume I > 1 and m > 1. In order to compute Sn{l,m), it is 
enough to consider only the following case: 

Wl W2 W3 '''27 

C(s, w) (f^) (/, m) = QQ^Tp^~PQQ^~Q ■ ■ ■ PP^~P 

where w — (wi, W2, ■ ■ ■ , ^27) G Z+' with wi, W2, ■ ■ ■ , ^27 > 1 and 7 > 1. For example, QPQP 
case is twi = ^2 = tfa = ^4 = 1 and 7 = 2. Then we have 

Z = W2 + W4 H h W27 

m = Wl + H h ^27-1 

As in the case of part (c), we obtain 
When Z A m = 0, this case is 

Sn{l,m) = 

It should be noted that 

r„(/,m) = s„(/,m) 

for any / and m. When I Am > 1, combining the above results (a) - (d) gives the proof of 
part (i). Proofs of parts (ii) and (iii) are trivial. So the proof of Theorem 3 is complete. 



Appendix B: Computation of the mth Moment 
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We begin with 



E 



J-\h/2 7r(l - x^)\/l 



2(3-2n)/2 

, dx — X lori 

/V2 7r(l-x2)Vl-2x2 TT 



where Z'^ is a weak hmit distribution of the rescaled Hadamard walk starting from 99 and 



sin" e 



It should be noted that 



(2 - y2)^i _ y2 y Jo l + cos2^ 
^^72 {2- (l + cos2^)}sin2"^ 



de 



'2n+2 



/■71 



2/. 



2n 



1 + C0S2^ 

(2n-l)(2n-3)---l tt 
2n(2n - 2) • • - 2 2 



Let J„ = -^2n/2". Then we see that for n > 



Jn+l Jji — 



23n+2 



f:) 



where 



Jo 



TT 



2V2' 



So we have 



-^2n — 2"'Jn — 2"" 



2V^ 



, 2-^ 
4^2 

"-^ 1 f2k^ 

k ; 



S 23fe+2 I 



Therefore for n > 1, 



k=0 



n-1 



TT 



1 /2A; 



In particular, 



2-^/2 
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